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Получены кинетические и полуэмпирические выражения для функции ползучести в физических уравнениях 
теории старения и кинетической теории прочности При сложном напряженно-деформированном СОСТОЯНИИ, 
а также для операторного модуля вязкоупругости и операторного модуля вязкого сдвига. 
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Введение. Современные несущие конструкции автомобильной, авиационной, ракетно- 
космической и других видов техники изготавливаются из композиционных материалов, позво- 
ляющих добиться реализации сложных конструкторских задач для достижения многофункцио- 
нальности, оптимизации массово-геометрических, экономических и других характеристик. При 
этом основной задачей проектирования является расчет длительной прочности несущих конст- 
рукций в условиях продолжительной эксплуатации. Трудность решения этой проблемы состоит в 
необходимости учета ползучести полимерных композитных материалов при сложном напряженно- 
деформированном состоянии в условиях многофакторного воздействия нагрузки. 

В связи со сложными условиями работы полимерных материалов в конструкциях установ- 
ление связи между напряжениями и деформациями при произвольном характере нагружения вяз- 
коупругого тела представляет собой одну из основных проблем механики полимеров [1-8]. Систе- 
ма уравнений, выражающая отмеченную связь, представляет математическую модель механиче- 
ских свойств материала. 

Для высоконаполненных полимеров такие модели прошли путь от дискретных до нели- 
нейных вязкоупругих моделей наследственного типа. Однако обобщение и анализ эксперимен- 
тальных данных, полученных для различных режимов испытаний материалов [1, 9, 10], показали, 
что принцип суперпозиции для высоконаполненных полимеров не выполняется, а построенные на 
его основе наследственные теории вязкоупругости не вполне пригодны для описания реального 
поведения материала. Кроме того, как правило, в процессе нагружения наполненного полимера 
изменяются его механические свойства. При повторном нагружении поведение материала зависит 
от характера предшествующего нагружения. Экспериментально установлено, что вязкоупругое 
поведение полимерного композитного материала при нагружении обусловлено вязкоупругостью 
полимерной матрицы и нарушением сплошности системы [1, 9, 10, 14]. В настоящее время не су- 
ществует методов непосредственного измерения степени накопленной внутренней поврежденно- 
сти рассматриваемых материалов [1]. 

Процессы старения полимеров и накопление повреждений в конструкциях при длительной 
эксплуатации дают основания применять модели и методы механики разрушений для исследова- 
ния этих процессов [15-17]. Однако имеющиеся методы и подходы не учитывают такие важные 
особенности развития поврежденности в полимерах, как диффузное проникновение легких фрак- 
ций через свободные поверхности. Технология полимеризации не исключает появления внутрен- 
них полостей или свободных поверхностей внутри материала. Наиболее активные изменения 
свойств, происходящие на этих поверхностях, в том числе и во внутренних микрополостях мате- 
риала, не учитываются в известных моделях полимеров. Для создания адекватных моделей на- 
пряженно-деформированного состояния конструкций из полимерных композитов необходимо учи- 
тывать поведение отдельных слоев вязкоупругопластических материалов, с учетом особенностей 
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временных изменений свойств этих материалов. Это обстоятельство требует учета кинетики хи- 
мических процессов в полимерах [12]. Множество факторов, характеризующих воздействие среды 
на состояние сложных конструкций, требует применения комплексных и обобщенных методов 
расчета, согласующихся между собой и с данными экспериментальных исследований. Созданию 
таких методов посвящено сравнительно небольшое число работ [18-20]. Проблема многофактор- 
ного анализа до сих пор остается открытой по причине лишь частичного рассмотрения вопросов, 
связанных с построением моделей комплексной оценки сложных систем. Полная модель прочно- 
стной надежности включает в себя, по меньшей мере, четыре составляющих, характеризующих 
свойства материала, конструкции, нагрузки и критерии прочности. Рассмотрению особенностей 
многослойных конструкций посвящены работы [4, 7, 11, 13, 19-23]. Однако практических алго- 
ритмов реализации методов расчета многослойных анизотропных конструкций намного меньше 
[11, 13]. Сложность создания вычислительных алгоритмов состоит в накоплении систематических 
и случайных ошибок при больших объемах вычислительных операций, отсутствии достоверных 
сопутствующих корректирующих и управляющих моделей, способных устранять погрешности вы- 
числений. Отдельную проблему составляет снижение временных затрат и оптимизация процессов 
вычислений. Кроме того, решение комплексных задач требует применения и комплексных мето- 
дов их решения, а также правильного определения условий взаимодействия этих методов, исходя 
из достоверной физической и математической интерпретации условий краевых задач. 

Таким образом, задача практической реализации расчетных методов исследования много- 
слойных композитных конструкций, учитывающих все рассматриваемые составляющие комплекс- 
ной модели прочностной надежности при сложном напряженно-деформированном состоянии, ос- 
тается в настоящее время не решенной. 

С этой целью рассмотрим способ перехода от модели ползучести материала при сложном 

напряженно-деформированном состоянии к модели вязкоупругости на основе принципа соответ- 
ствия путем определения выражений для операторного модуля вязкоупругости и операторного 
модуля вязкого сдвига. Для получения указанных выражений определим функцию ползучести с 
помощью кинетических уравнений [12] и эмпирических констант [10]. 
Уравнения ползучести на основе теории старения и кинетической теории при сложном 
напряженно-деформированном состоянии. Полная система уравнений теории ползучести, 
так же, как и в теории упругости и пластичности, включает три основных группы уравнений: рав- 
новесия или движения, уравнения физические и геометрические. Рассмотрим уравнения ползуче- 
сти на основе теории старения и кинетической теории для полимерных композитных материалов 
(ПКМ), например для конструкций летательных аппаратов, находящихся в условиях длительной 
эксплуатации и стационарного термосилового нагружения (ТСН). 

Уравнения равновесия в тензорном виде: 


до’, 
>А+Х, =0, 1,1 =1,2,3. (1) 
дх, 
Геометрические уравнения Коши: 
1 ди, ди, 
== —“+—^ |, 5/1 =Ь2,3. 2 
О: 0 ы (2) 


Здесь о,, &,, и, - соответственно компоненты тензоров напряжений, деформаций, пере- 


мещений; х, - переменные осей координат; Х, — компоненты объемных сил в направлениях х,. 


Физические соотношения в предположении, что мгновенные деформации являются упру- 
гими, кривые ползучести — подобными, что существует степенная зависимость деформаций пол- 
зучести от напряжений, удобно записать в следующем виде: 


1 в 
21 == 6 —н(с,› +03 )] +8; о ‘(0 (о 56») 
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1 Зла 

&> и -и (Оз +6445; ‘(0 (°..-9.,}; 
1 Е 

533 = -и (о +0 )| +: ый '9(0 (с; -0.); 


б ы , 
&5 = о 30" 'О(/а.›; 


с ыы 
= = +30" 'О (дб; 


вз= 9 +397 Ось. (3) 


1 
Здесь =, - температурные деформации; с., то +0, +0,,) - среднее напряжение; 


Е 

п 

7 
маций ползучести. 

Функцию ползучести и температурные деформации определим с помощью кинетических 

уравнений [12]. В соответствии с [10], деформация ползучести может быть определена в виде 


произведений функций: 


О(В = — функция ползучести; с, — интенсивность напряжений; 5” — интенсивность дефор- 





5" (©, =у(6)0(0. (4) 
В практических расчетах функцию \у(с) считают степенной: 
(с) =0", (5) 


где и>1 — постоянная для данной температуры и материала величина. 
Используя полученные выше выражения для деформаций ползучести, определяем выра- 
жения для функций \у(о),О(г). Экспоненциальную функцию =” можем записать через эмпириче- 


ские константы [10]: 


1 
&=—АВое-“ = Воем =\(0о(В. (6) 
Е Ц 
Тогда, в соответствии с выражениями для напряжений, получим: 
(с) =Во. ; (7) 
1 1 
ОА "=". (8) 
Е Ц 


В кинетическом виде искомые функции можно определить с помощью уравнений [12]: 
1 ГА 
в" =-ош— =\(ДО(. (9) 
Ц бо 
Сравнивая полученные уравнения (8) и (9), запишем аналогичные по виду кинетические 


функции: 
(с) =с; (10) 


О-о. (11) 
т 


0 
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В соответствии с теорией старения, сомножитель с” "'©(Р) в выражении (3) определяется 


соотношением: 
П 
в 5. 
(==. 
[©] 


7 





Тогда согласно (6)-(11) сомножитель (12) можно записать: 
— через эмпирические коэффициенты: 


П 
- Е 6 ы 9%, 
оО =-=—Ае” = Ри ты 
о, Еб; по; 
— согласно кинетической теории: 


оо =А шв. 
т 5 


(12) 


(13) 


(14) 


Используя полученные выражения (13), (14), а также кинетическое выражение для тем- 
пературных деформаций [12], можем записать для каждого фиксированного момента времени # 


систему физических уравнений теории старения в полуэмпирическом виде: 


1 3 Ве, 
= = [бы = (©, +5; )|+аАТ+ 2” м (с -в„)}; 


1 3.Вбо 

>> = 6 -н (о, она бе " (5» 6 
1 ЗВо. , _ 

& == 6» и (< ОА ь (.; -с.,); 


1 
вв =6ь | + Зем ва . 
То; 


а - +34" в Е 
С по; 





1 
вв =6в | 2+ Зем в р 
Тс; 


или, согласно выражению (14), полученному на основе кинетического подхода 


1 [от 3 1 
и =—|6и-и(о»›+6.. )|+ш— | —— +9, -с©.,)- |; 
и — и-Н(°>> зз ы ( И „= 





ЕС 2 
1 ВЕ 1 
&> == 6» —=н (с. а [с о =.) р 
1 (| одт 3 1 
&зз 6» = (оп УБЛ: | с 0195 -«„)= | 
вые (чзнь: 
с то 
1 1 1 
#3 =6:| —+3—Ш— |; 
23 (5 а я 


(15) 


(16) 
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Анализ полученных уравнений показывает, что в физических уравнениях для деформаций 
сдвига последний сомножитель играет роль операторного модуля сдвига С’ теории вязкоупруго- 
сти, соответствующего модулю сдвига С’ в теории упругости, по аналогии с определенным ранее 
операторным модулем упругости Е. 

В соответствии с полученными кинетическими и полуэмпирическими уравнениями (15) и 
(16) теории старения, можем записать выражения для операторного модуля сдвига: 

— в полуэмпирическом виде: 





—1 
= 1 с 
2 [= зе" Вбь |, (17) 
а тб; 
— на основе кинетического подхода: 
—1 
р 1 1 ГА 
= [ зи о (18) 
а то 
Перепишем кинетическое уравнение для одноосного деформирования в раскрытом виде: 
1 Е ОТ ГА ГА 
и =—0, + я +9, —Шш—. (19) 
и бо т 6 


После преобразования выражения (19) по аналогии с выражением для напряжений кине- 


тической теории [12] получим: 
Е [аАТ р 
$1 =50, + —| —— +06, — |. (20) 
п С Ц 
Учитывая, что второе слагаемое в скобках выражения (20) представляет собой напряже- 
ния, обусловленные тепловой составляющей, запишем: 


= = +в (Г) +6 в. (21) 
и т 5 

Тогда для обобщенных напряжений, обусловленных действием тепловой и механической 
составляющей ТСН, получим обобщенное уравнение связи деформаций и напряжений при одно- 
осном деформировании: 





1 
ев (22) 
Е о 
или после преобразования: 
01 11 
& =—(1+-—ш—#). (23) 
. Е т 5 
Уравнение для величины, обратной операторному модулю упругости, примет обобщенный 
вид: 
1 1 1 
В (24) 


Е; Е, Ц бо 
Таким образом, на основании принципа соответствия теории упругости и вязкоупругости 
справедливы физические соотношения между напряжениями и деформациями ползучести для 
реологических материалов, выраженные через операторные модули вязкоупругости Е и сдвига 
С, ‚ полученные в кинетической и полуэмпирической форме (17), (18), (24): 








57а =9#, 1=1,2,3, }=1,2,3 2 
== и би и,2, =Ь4л,2. (25) 
Е: С, 
Полученные соотношения (11)-(25) позволяют проводить анализ ползучести полимерных 
конструкционных композитов при сложном напряженно-деформированном состоянии с помощью 
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методов теории упругости, при использовании полученных выражений для операторных модулей 


вязкоупругости и вязкого сдвига. 
Общий алгоритм анализа ползучести полимерных композитов при сложном напряженно- 


деформированном состоянии представлен с помощью блок-схемы (рис.1). 
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Рис.1. Блок-схема алгоритма анализа ползучести полимерных композитов 
при сложном напряженно-деформированном состоянии 


В соответствии с этим алгоритмом на первом шаге составляется исходная система уравне- 
ний упругости для сложного напряженно-деформированного состояния, в которую входят урав- 
нения вида (1), (2), (25). Далее, на втором шаге, в соответствии с методами [11], записывается 
разрешающая система уравнений относительно функций распределения напряжений по толщине, 
в которых выражаются условия контакта между слоями. На 3-м и 4-м шагах осуществляется ввод 
исходных данных о Ффизико-механических (вязкоупругих), теплофизических и массово- 
геометрических характеристиках конструкции, а также исходные данные о нагрузках, выражен- 
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ные через функции напряжений, перемещений и температуры. Указанные исходные данные ис- 
пользуются для разложения функций параметров нагружения и напряженно-деформированного 
состояния в двойные тригонометрические ряды по поверхностным координатам. На шаге 5-м оп- 
ределяются коэффициенты рядов Фурье. На шаге 6-м вычисляются гармонические компоненты 
функций разложения для Фурье-анализа. Одновременно на шаге 7-м осуществляется разложение 
производных по времени в конечные разности. Полученные результаты разложения функций ис- 
пользуются для разделения переменных в разрешающей системе уравнений, полученных на вто- 
ром шаге алгоритма. В результате разделения переменных и вычисления производных от простых 
тригонометрических функций, выполненных на шаге 8-м, приходим к системе обыкновенных 
дифференциальных уравнений на шаге 9-м, для решения которой, на шаге 10-м, задаются пара- 
метры расчета, определяющие число гармоник по поверхностным координатам, число точек ин- 
тегрирования по толщине и дискретной ортогонализации. 

Кроме этого, на шаге 11-м, осуществляется ввод данных о начальных и граничных услови- 
ях для решения краевой задачи. В результате решения задачи, на шаге 12-м, определяются иско- 
мые полные функции напряжений и перемещений путем двойного суммирования волновых функ- 
ций в соответствии с методом анализа Фурье. Решение осуществляется до полного окончания 
суммирования по всем волновым числам. На этом алгоритм завершается (шаг 13), и результаты 
решения выводятся на экран монитора (шаг 14). 

Проверка достоверности используемой модели показала точную сходимость результатов 
расчета тестового варианта для изотропной тонкостенной оболочки (рис.2, а) с результатами того 
же варианта по классической теории с «наложением» решения по трехмерной теории в области 
краевого эффекта [11] (рис.2, 6). 


О’, Н/м?.10-7 О’,Н/м?:107 О’, Н/м?'1077 





в) 

Рис.2. Распределение окружных и меридиональных напряжений по наружной и внутренней поверхности оболочечной 
конструкции: а - из металла АМг-6, полученное с использованием классической модели [13]; 6, в - из металла АМг-6 и 
композитного материала СВМ + ЭДТ-10, соответственно, полученное с помощью представленной трехмерной модели 
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Расчет проводился для металлической оболочки из АМг-6 с размерами: [/А=0,5, 
/К=0,25, фо=0,25Г, (рад), где Г, К — длина и радиус оболочки; (1, фо — линейный и угловой разме- 
ры локальной области. Результаты расчета окружных и меридиональных напряжений наружной и 
внутренней поверхностей оболочки представлены на рис.2, а, 6. 

С помощью предложенной модели были проведены исследования многослойных анизо- 
тропных оболочек из композитного материала СВМ + ЭДТ-10. 

Результаты расчета напряженно-деформированного состояния (НДС) при нагреве локаль- 
ной области для рассмотренного выше варианта нагружения композитной оболочки приведены на 
рис.2, в. 

Сравнение полученных результатов расчета показывает общее сходство картины распре- 
деления напряжений в композитных и металлических оболочках, однако, напряжения в компо- 
зитной конструкции в 3-3,7 раза меньше, чем в металлической. В этом проявляется преимущество 
композитов перед металлами при использовании их в качестве конструкционных материалов в 
различных областях техники. 

Заключение. Приведенная математическая модель, а также методы и алгоритмы ее применения 
позволяют прогнозировать кратковременную и длительную прочность современных и перспек- 
тивных несущих конструкций из полимерных композитных материалов при многофакторном воз- 
действии нагрузок и условий эксплуатации на основе численных методов расчета напряженно- 
деформированного состояния многослойных оболочек с использованием анализа Фурье. 
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